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Abstract 
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I. INTRODUCTION 



Few-body systems with resonant interactions show universal properties independent of 
the details of the interaction at short distances [1]. If the scattering length a is much larger 
than the range of the interaction tq, the underlying interaction can be approximated by an 
effective theory with contact interactions. In the two-body system, the universal properties 
include the existence of a shallow dimer with binding energy E2 ~ —h'^/{ma^), if a is 
positive. In the three-body system, the Efimov effect then generates a geometric spectrum 
of three-body bound states [2]: 

E^^^ = -{e-'^/^^rh'Kl/m, n = 0,1,2,... (1) 

where sq = 1.00624 . . . and is the binding momentum of the lowest Efimov state. This 

spectrum shows a discrete scaling symmetry which corresponds to a renormalization group 
limit cycle [3]. Recent theoretical studies have shown that these universal properties also 
hold for the four-body system. In particular, there is a pair of universal four-body states 
associated with every Efimov state [4, 5]. Moreover, universal cluster states have been 
predicted in five-, six- and possibly higher-body systems as well [6-9]. Using Feshbach 
resonances, the two-, three-, four- and five-body states have been confirmed in a ultracold 
atomic gases in a variety of experiments with different atom species [10-12]. 

These experiments were carried out in a regime where the influence of the trap on the few- 
body spectra could be neglected. However, the trap also offers new possibilities to modify 
the properties of few-body systems. In particular, a narrow harmonic confinement modifies 
the spectrum. Microtraps with a small number of atoms offer the opportunity to study the 
transition from few- to many-body systems in a controlled environment. 

Many studies have focused on fermions in the unitary regime in a harmonic trap where the 
Efimov effect does not occur (see Ref. [13] for a review of this work). Here, wc focus on the 
case of spinless bosons. The 2-body problem of spinless bosons with contact interactions in 
an isotropic harmonic trap was solved analytically by Busch et al. [14]. The corresponding 
3-body problem was first solved by Jonsell et al. [15]. Werner and Castin calculated the 
complete 3-body spectrum in the unitary limit of infinite scattering length and provided a 
semi-analytic solution [16]. More recent work has focused on the universality of the two- 
body spectrum [17, 18], the extension to heteronuclear systems [19] and anisotropic traps 
[20, 21], and the interpretation of the scattering- length dependence of the level structure in 
terms of the Zel'dovich effect [22] . 

In this work, we investigate universal few-body physics in harmonic confinement for up 
to six bosons using an effective theory framework. Our strategy follows Stetcu et al. [23- 
25] , where an effective theory for short-range nuclear forces in the framework of the no-core 
shell model was formulated. The effective interactions were defined within a finite model 
space with a cutoff N on the basis functions. In previous work [26], we have performed a 
systematic study of the spectra of three- and four-body systems in the vicinity of the unitary 
limit at leading order in the effective theory. In this paper, we extend our calculations to six 
bodies. The running of the coupling constants for two- and three-body interactions with the 
cutoff N is investigated in detail. Moreover, we introduce smeared contact interactions. We 
compare the convergence of the effective theory as the size of the model space is increased 
for exact and smeared contact interactions and discuss different strategies to remove the 
regulator. Finally, spectra for A-boson systems with up to six bodies are calculated using a 
combined extrapolation in N and the smearing parameter e. 
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II. THEORETICAL FRAMEWORK 



We investigate few-boson systems with resonant two-body interactions characterized by 
a large S-wave scattering length a in a harmonic trap using an effective theory with contact 
interactions. The contact interactions lead to ultraviolet divergences which have to be 
renormalized by fixing one two-body energy (or, equivalently the scattering length a) and 
a three-body energy [3, 27]. Since the trap only modifies the infrared properties of the 
system, the renormalization of the contact interactions is the same as in free space. We 
stay at leading order in the effective theory. (See Refs. [28, 29] for an explicit calculation 
of effective range corrections for two-component fermions.) Higher-order corrections are 
suppressed by kR where k = a/ m\E\ is a typical momentum and R is the range of the 
interaction. For broad Feshbach resonances R is given by the van der Waals length /vdw- 
For trapped systems, there are also corrections of order -R//5, where /3 = ^/Uj (mu) is the 
oscillator length. As we will discuss below, there are additional corrections if the two- and 
three-body renormalization energies E2 and are very different. 

We implement this effective theory in a Hamiltonian framework. For a system of A bosons 
in an isotropic harmonic oscillator potential (HOP) at leading order the Hamiltonian can 
be written as 



^ = E fl!^ + 'o^^m'] + 9''' E + 9''' E w,, , (2) 



j=l \ / i<j i<j<k 



where V and W are two- and three-body contact interactions, m is the mass of the bosons, 
and u is the trapping frequency. In the following, all lengths are rescalcd with the oscillator 
length /3 = ^JhjimuS) for convenience. Since the interactions are Galilei invariant and 
the HOP allows to separate the center of mass motion, the problem is transformed to an 
(A — l)-body problem in Jacobi coordinates Sn and a one body problem in R^ where 



1 ^ . ^ 1 



-,•=1 pVA 



with the property Y^f^^ = (Zl^=i l^'np + 

The Hamiltonian (2) has to be regularized since it is not self-adjoint. Hence the Hilbert 
space is restricted to the model space consisting of the linear hull of oscillator states 
{(S)t=i\nilimi)} with eigenvalues j2t=i ^C^^i + h + < huj{N + {A - 1)|) [30]. In 
order to match the results in the model space to the full Hilbert space results, the coupling 
constants g^'^^ and g^^^ are renormalized (see Sec. II B). We identify with the ultraviolet 
cutoff of our effective theory. The couphng constants g^"^^ and g'^^^ are fixed by matching to 
a two- and a three-body energy, respectively, and run with the cutoff A^. In free space one 
usually uses a momentum cutoff A to regularize the resulting integral equation and finds a 
linear dependence of observables on 1/A in first order of perturbation theory. Because of the 
relation p oc \/E between momentum p and energy E, we expect errors in the eigenenergies 
in the trap of order 1/ ^/WTJA^^Tj3/2. 

In contrast to the momentum cutoff A used in free space, the regulator A^ also implies an 
infrared cutoff. In the literature, there are two definitions of the IR-cutoff for an oscillator 
basis. Both the expression f3/^N+{A - 1)3/2 [31, 32], which quantifies the maximum size 
of structures that can be captured in the given basis, or the oscillator length /3 itself [33] have 
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been interpreted as an IR-cutoff. In our calculations, the harmonic oscillator represents a 
physical trap and thus /3 is a physical parameter. The trap acts like a finite box that confines 
the system. In this case, only the first definition is sensible and we expect errors of order 
1 / -^/N -\- {A — 1)3/2 from the infrared cutoff. These errors have the same scaling behavior 
with N as the errors from the ultraviolet cutoff and vanish for N ^ oo. 

In our previous work, we found that the convergence to the results in the full Hilbert 
space is slow due to the irregularity of the ^-distribution in the parametrization of the 
genuine contact interactions [26]. For the cutoff dependence to become insignificant, values 
for 1/y/N of order 0.1 or smaller have to be reached. For four- and higher-body states such 
low values can in general not be obtained. The eigenenergies are estimated by extrapolation 
of results for finite N to N ^ oo. 

Here, we investigate a second possibility. The convergence can be improved by smearing 
the contact interactions. We keep the separability of the interaction and approximate the 
contact interactions by narrow Gaussians with width e. More specifically, the separable 
two-body contact interaction 

{s\V\s') = 5^'\s)6^'\s- s') = S^'\s)6^^\s') (4) 

is substituted by 

{s\v4s') = (27re^)-^/^e-l^1'/(2^') (27re^)-^/2e-l^""l'/(2^') , (5) 
and the three-body contact interaction 

{si, s-,\W\s',, 4> = 5('\s-l)S^'\s[)S^'\s-,)5^'\s',) (6) 

is replaced by 

{si, s-'2\W,\s[, 4) = (27re2)-6 e-(l^^^l'+l^^^l'+l^""il'+l^"*^l')/(2^') . (7) 
The case of contact interactions is then recovered in the hmit e — > 0. 

A. Matrix elements 

The eigenenergies are determined by diagonalization of the Hamiltonian matrix in an 
angular-momentum-coupled basis of oscillator states. For contact interactions as well as 
smeared contact interactions one can easily find the matrix elements for the oscillator basis 
functions 

0nihmi(Sl) = RnJi{\Si\)Yi^rni{^) = (si|ni/imi> . (8) 

For two-body contact interactions, we have 

{nihmi\V\n[l[m[) = <^„io(O)0„;o(O) Si^qSi'^q , (9) 



where 



WO)- ' (J>+Jl_V''_({2n + mY' .^Q. 
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For the smeared two-body contact interaction, we obtain according to Appendix B 

1/1 — e^\"i+"'i 

{nihmi\V,\n[l[m[) = ^-^ ^ ^2^3 ( 1 + g2 j 4>mo{0)4>n{o{0)6i,oSi'^o . (11) 

The three-body interactions are determined in a coupled scheme, i.e. both angular mo- 
menta are coupled to a total angular momentum with Clebsch-Gordan coefficients, 

[(l>nih{Si) ® 0„2/2(^2)]ij = {si,S2\nih,n2l2,LM) . (12) 

.-2 = ^(^1+X2-2X3), 

^ -^;^(^i + ^2 + ^3) • 

Figure 1. Set of Jacobi coordinates for the three-body system. 

We find for the three-body contact interactions in Jacobi coordinates Si and S2 
{nih,n2l2,LM\w\n[l[,n'2l'2,L'M') = 0„io(O)0„;o(O)0„2o(O)0^o(O) V^r.o^i^o^i^o , (13) 
and finally for the smeared three-body contact interaction 

(711/1,^2/2, LM\W,\n[l[,n2l2, L'M') = 

(1 + e^)6 V I + 0nio(O)0n2o(O)0„'^o(O)0„^o(O)5«io5/io5«^o5i2O • (14) 

B. Renormalization 

Next we discuss the renormalization of the Hamiltonian, Eq. (2). For convenience, all 
energies are given in units of hij. The coupling constants g^'^^ and g^^^ depend on the 
regulator A^. At each value of A^, i.e. in each model space, their values have to be fixed by 
one two-body and one three-body observable. All other observables are then predictions. 

As the two-body input parameter, we use an eigenenergy E^"^^ in the two-body spectrum 
with L = 0. Using the analytical solution derived by Busch and collaborators [14], 

(3 ^^ r(3/4-E(^)/2) 

a ""^ r(l/4-E(2)/2) ' ^ ^ 

this energy E^'^^ is directly related to the S-wave scattering length a. 

The relation between g^'^'^ and E^'^'' for a given A^ can be found analytically by exploiting 
of the separability of the interaction. We define 

{</)„o(0) , for contact interactions, 

(1 + e^)^^/^ ^Y^^ 0no(O) , for smeared contact interactions. ^ ^ 
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As discussed in Sec. II A, the matrix elements in the 2-body sector for L = Z = are 
given by 

{n \H('^ I n') = (^2n + ^ Sn,n' + ■ (17) 

In the model space corresponding to N, the solution \1pEC2}) of the eigenvalue problem with 
the eigenvalue is expanded in oscillator functions 

N/2 

|^^(^)>=E^?1^>- (18) 

n=0 

Applying the Hamiltonian to Eq. (18) and projecting on the oscillator state \k), we obtain 

2^ + I) of + E 9^V.Jn.cr - Ei%r ■ (19) 

^ n=0 

Solving for cf^^' and reinserting the result in Eq. (19), we find the running of the coupling 
constant g^'^^ with N: 



The three-body coupling constant g^^'' is fixed by an energy of the three-body system E^'^\ 
The solution \aN) of the three-body system without three-body interaction for a given value 
of N is: 

1"^)= E ZNiaN;nili,n2l2,L)\nili,n2l2,L) , (21) 

nih,n2l2,L 

with energy eigenvalue D]\f{aN)- Note that the eigenstates are degenerate in the total 
angular momentum projection M. With these states, we calculate the matrix elements for 
the complete Hamiltonian including the three-body interaction 

+ g^^\N) ( E Zn («iv; niO, n^O, 0) ( E ^n (a^; ^'lO, n',0, 0) /^^ . (22) 



ni,n2 



Requiring the energy of the state jaAr) to be E^^^ for all A^, the renormalization condition 
for the three-body coupling follows: 



^(3) (TV) DN{aN)-E(^) 



Note that the expansion coefficients Z^v (ccjv; riiO, ^26, 0) as well as the eigenvalue Dj^{aN) 
explicitly depend on A^. The full spectrum in a model space for given A^ can then be 
determined with these coupling constants by diagonalization of the Hamiltonian matrix. 
Our general strategy for calculating the Hamiltonian matrix for a system of A identical 
bosons is described in Appendix A. 
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C. Running of g^^^ and g^^^ 



In this subsection, we study the running of the couphng constants g^"^^ and g^^^ for contact 
interactions in detail. Stetcu et al. [28] rewrote the sum in Eq. (20) in terms of F-functions 
and the generahzed hypergeometric function 3F2. They found an exphcit relation for g^'^\N) 
from which the behavior for large can be obtained. 

Here, we provide an alternative derivation of the behavior of g^'^\N) for large values of 
N using an integral representation of the sum in Eq. (20). For E^^^ < 3/2 each term in 
the sum is positive and the denominator grows monotonously. We examine the behavior 
for very large N. With the Euler-Maclaurin formula the sum can be approximated as an 
integral 

V <PloiO) _ 1 F(a; + 3/2)da; 

^2n + |-£;(2) TT^^J, F(3/2)F(a; + l)(2a;+|-E(2)) ^ ' ^ ^ 



^(2) (AT) 



where (t)no{0) was substituted with Eq. (10). The quotient of Gamma functions can be 
expanded in x, 

r(. + 3/2) _ 1 ^^o(l). P5) 



F(3/2)F(x + l) F(3/2)^ 
Inserting this expansion in Eq. ((24)) and integrating, we find 



Thus the coupling constant vanishes as l/y/N. Identifying vA^ with the momentum cutoff 
A, this is consistent with the renormalization in free space [3]. We thus expect the leading 
errors from finite N in our effective theory to scale with 1 / \/N . 

For E^"^^ < 3/2, the coupling rapidly approaches zero as is increased. In the case that 
E^'^^ > 3/2, the terms in the sum in Eq. (24) are negative at first until n > (E^^^ /2-3/4). The 
coupling g^'^\N) as a function of iV G M thus develops a minimum for AT (£;(2)/2-3/4) 
and has a pole as A^ is increased further. For even larger A^ it approaches zero as well. 

For smeared contact interactions an additional damping factor appears. For Eq. (24), we 
have 



E 



^(2) {N) (1 + £2)3 ^ (2n + I - Ei^)) \l + e\ 

^ L^r/^ r(x + 3/2)d. (l^X\oil) (27) 

7^3/2(1 + ^2)3 70 r(3/2)r(x + i)(2x + |-£;(2)) Vi + eV ^ ^'-^ ^ 

The integral now converges to a constant ior N ^ 00 exponentially at a fixed value of e. 
Thus, the coupling constant converges to a finite number. Note, that the larger e is, the 

faster the coupling constant converges and the energy spectrum becomes independent of A^. 
This result reflects the additional regularization of the contact interaction by the smearing. 

In the inset of Fig. 2, we illustrate the behavior of the running coupling constant \fNg^'^^ 
for e = 0. We show \fN g^'^^ in the unitary limit for the renormalization energies £'^2) — 0.5 
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(ground state) and E^'^^ = 16.5 (8th excited state). In the case of E^'^^ = 0.5, the running 
couphng g^'^^ has aheady converged to a constant value for small values of and shows 
no structure. For E^'^^ = 16.5, the situation is as described above. The coupling constant 
approaches zero for around 14 where the denominator of the right hand side of Eq. (24) 
changes sign. Moreover, it changes from large positive to large negative values around 
N = 21. This behavior can be understood by looking at the spectrum for E^"^^ = 16.5 
shown in Fig. 2. For small values of N, the model space is not large enough to describe 
the ground state adequately. The behavior of the coupling constant a/ZV is exactly such 
that this new state enters the renormalized spectrum from minus infinity keeping the other 
states unchanged. In the continuum case, a similar behavior is observed for the three-body 
spectrum [3, 34]. For larger values of A^, when the model space is large enough to describe 
the ground state, the running coupling g^'^^ approaches the same value for both cases. 




1/(N+3/2)^'^ 



Figure 2. Energy spectrum in the two-body sector for E^"^^ = 16.5. The dots on the left show the 
exact energies. The inset shows the running coupling constants \/iV g^"^^ in the unitary limit for 
renormalization energies = 0.5 and £^(2) = 15,5^ 

For the coupling constant g^^^ the situation is more complicated. The eigenvalues D(a) 
as well as the eigenstates now depend on the cutoff parameter A^. Therefore it is not 
straightforward to derive the leading-order behavior of g^^^ analytically. Werner et al. have 
provided a semi-analytic solution for the three-body problem in a harmonic trap in the 
unitary limit [16]. We use their results for the energy spectrum to benchmark our three- 
body calculations. 

In the left panel of Fig. 3, we show the spectrum in the unitary limit for a three-body 
renormalization energy ^ 1.76 as a function of A^. At A^ = 16 a new three-body 
state enters the spectrum of the model space from minus infinity. For large A^ this state 
approaches the exact eigenenergy at —4. In the inset, the corresponding coupling constant 
g'^^^ is shown. As in the two-body case, the coupling constant diverges, here around A^ = 16, 
and changes from large positive to large negative values. The three-body spectrum and the 
coupling constant in the continuum case show the same behavior [3]. The phenomena are 
the same as in the two-body case discussed above. In the right panel of Fig. 3, we present 
the energy spectra for the renormalization E^'^'^ ^ 16.14. This renormalization belongs to 
the energy spectrum in the full Hilbert space: ■ ■ ■ , -5319.28, -10.29, 1.4, 3.68, 5.82, 7.92, 
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Figure 3. Left panel: Energy spectrum of Efimov-like states in the three-boson sector in the unitary 
hmit (E'^^^ = 0.5) for E^^') w 1.76. The dashed hnes indicate the exact energies [16]. Right paneh 
Energy spectrum for Efimov-hke states in the unitary hmit E^'^^ = 0.5 with E^'^'^ ~ 16.14 for contact 
interactions. The running couphng constants N'^ g^^^ corresponding to both cases are shown in the 
insets. Poles are indicated by vertical solid lines. 



9.99, 12.05, 14.10, 16.14, ■ ■ ■ . The states with energies lower than 1.4 cannot be found inside 
the feasible model spaces for this renormalization energy. Remarkably, the states approach 
the exact energies but at specific values of the cutoff marked by vertical lines the spectrum 
rearranges. The smallest value of N corresponds to the appearance of a ground state at 
energy 1.4. 

At the two largest cutoff values marked, there are small discontinuities in the energy. 
These discontinuities are artifacts of the renormalization method and do not correspond 
to new states entering the spectrum. Since we work in a finite model space, not only the 
Efimov-like states but also the universal states depend very weakly on the three-body inter- 
action. In the special case when for a given cutoff N the renormalization energy E^'^^ exactly 
coincides with the eigenvalue of a universal state corresponding to another spectrum with 
renormalization energy our computer code erroneously renormalizes to the spectrum 

characterized by E^^\ The energy values in the neighborhood of such points should therefore 
be discarded. 

In the inset of the right panel of Fig. 3, the behavior of the corresponding coupling 
constant is shown. The coupling constant has a salient behavior at three positions indicated 
by the vertical lines. At the smallest value of N where the ground state enters into the 
spectrum from below, the coupling changes from large positive to large negative values. The 
approximate pole which is expected at this position due to the discrete nature of the cutoff 

is strongly distorted on the right wing due to finite cutoff effects. At the two larger 
values of A^, there are discontinuities due to the renormalization artifacts discussed above. 
Moreover, for even larger A^ the coupling approaches an approximate pole corresponding 
to the addition of a new ground state at energy —10.29. This pole is not reached in our 
calculation; it would require a larger A^. 
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D. Error Analysis 



There are various sources of errors in our calculation. In this subsection, we perform a 
detailed error analysis. 

First, there are corrections due to the ultraviolet cutoff parameter A^. Our considerations 
in subsection II C showed that the finite model space also implies an infrared cutoff which 
vanishes as N is increased. The errors due to both cutoffs show the same scaling behavior in 
A^. Therefore, we expect corrections in the energy eigenvalues of order 1/ ^^/N~^^JA'--T)S/2 
for large N. The shift of by {A — 1)3/2 under the square root takes into account the 
zero-point energy of the free A-hodj system. We can extract the scaling behavior of these 
corrections from the error analysis introduced by Lepage [35]. In Fig. 4, we show the 
deviation of the lowest three-body energy eigenvalues from the exact values in the unitary 
limit E^'^^ = 0.5 with three-body renormalization energy E^^^ = —1. The double-logarithmic 
plot shows a linear dependence of the energy differences on log(A^ + 3) for large values. 
From a linear fit to the five largest values of A^, we find a slope of s ~ —0.6 for the first 
four Efimov states above the state used for renormalization. The small difference to the 
expected value of s = —0.5 could be due to contamination from higher order corrections. 
As a consequence, the cutoff dependence for contact interactions is in agreement with power 
counting arguments based on identifying the continuum cutoff A with 1 / a/A^ + (A — 1)3/2. 
A similar power law dependence of the leading corrections to three-body energies for contact 
interactions was observed by Furnstahl et al. [36]. 




0.1 I ■ ' ■ ■ ■ — ' — ■ — ■ — I 

10 100 

N+3 

Figure 4. Corrections in energy eigenvalues as a function of A + 3 with A = 3, E^^^ = 0.5, 

In the case of smeared contact interactions, we find an exponential dependence on \/N. 
In the two-body sector, we consider model spaces with cutoffs up to A = 700. In a Lepage 
plot of log{\EN — Eoo\) we find for large cutoffs a slope of about —0.5. This implies that for 
large N the calculated energies behave as 

En^E^ + cie-^^^ . (28) 

Thus the smearing changes the leading corrections from an inverse power law to an expo- 
nential behavior. 
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The convergence properties of variational calculations using basis expansions have been 
studied already in the 1970's [37, 38]. More recently, the convergence properties of ab initio 
calculations of light nuclei in a harmonic oscillator basis were investigated [33, 39]. Due 
to their singularity, we expect the contact interactions to behave quite differently, but it is 
interesting to compare our results for smeared interactions to Refs. [33, 39]. However, one 
has to keep in mind that we consider a physical trapping potential instead of a mere basis 
expansion. Moreover, in our calculation the effective interaction is renormalized at each N 
in order to keep the energies of a given two- and three-body state fixed. 

Our results show the same exponential dependence on the ultraviolet cutoff and on the 
infrared cutoff A^c oc l/^jN + 3/2 as observed by Coon et al. for the Idaho N3L0 poten- 
tial [33]. Furnstahl et al. [39], in contrast, observed a Gaussian dependence on the UV 
cutoff 

E^^E^ + Aoe-^^(^)' , (29) 

using similarity renormalization group (SRG)-evolved chiral interactions, but their depen- 
dence on the infrared cutoff is consistent with ours. The difference might be due to the use 
of SRG-evolved interactions, but this question requires further study. 

If we want to extract the energies for genuine contact interactions from smeared contact 
interactions, there are also errors due to the width parameter e which generates a finite 
interaction range. The corresponding two-body problem in free space is separable and 
can be solved analytically. In our calculations, we determine the coupling constants from 
matching to energy levels in the oscillator. One two-body energy level E^'^^ is kept constant in 
each model space and the coupling constant g^'^^ is determined from this matching condition 
(cf. Eq. (20)). We find: 

N/2 ( l-e'V'^ (2n+l)!! 



^-^ ((. . .) V.) - 3)-s.. (I ! - I - Ifef ) ' 

where 2Fi(a, 6; c; z) denotes a hypergeometric function. The limit N ^ oo yields the coupling 
constant in the full Hilbert space. 

We use e as a smearing parameter and are interested in the behavior for small e. Ex- 
panding the first and second line of Eq. (30), we find: 

£\e) = -2.i . - ^ ' - 2.i (l + 4B» + ^^^^w) + ' '^'^ 

U 2 ^ ^ U 2 ^ 

g^ie) = co(7V, E^'^) + c^iN, E('^)e' + 0{e') . (32) 



Remarkably, the results for finite and infinite N differ fundamentally since the hmits e — > 
and A?" — > oo do not commute. For N ^ oo, the coupling constant depends linearly on 

e in leading order and vanishes for e — )■ 0. For finite N, the leading term is independent 
of e and the first correction is of order e^. The constant term in e is a consequence of 
regularization with N but the vanishing of the correction linear in e is unexpected. The 
expansion parameter co{N, E'^^)) vanishes like 1/VN for A?" — >■ oo. 
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The corresponding two-body energies show a similar behavior. First, the value of (7(2) 
is fixed for a given e and E^'^\ The other energies E^'^^ can then be calculated numerically 
using a root-finding algorithm. For N = 00 and small e the energies Ef^(e,E(2)) depend 
linearly on e. In order to extract the results for a zero-range contact interaction from a 
calculation with smeared interactions, a linear extrapolation in e is thus appropriate. For 
finite A^, however, the dependence of the excitation energies Ef'\e, E^'^\ N) on e is different. 
For e — 0, we find \Ef\e, E^'^\ N) — Ef\{], E^'^\ N)\ oc and an extrapolation in e does not 
reproduce the result for a zero-range contact interaction. The e-dependence of the energies 
is illustrated schematically in Fig. 5. In order to extract the contact interaction results from 
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Figure 5. Schematic dependence of energies on the smearing parameter e for finite N (dashed hne) 
and infinite N (sohd hne). 



smeared interactions, it is therefore important to extrapolate in first (or use numerically 
converged values) and then extrapolate to e = in a second step. 

The variation of the energies with e can be interpreted as a dependence on the effective 
range r^. of the two-body interaction. The effective range parameters in terms of the effective 
coupling constant g'^'^^ in the continuum case are (See Appendix C): 

1-^ + ^ and r-^-^ (33) 

These relations are valid in the full Hilbert space but not in the restricted model space 
characterized by finite A^. Analog to the case without a trap, our numerical results and 

(2) 

Eqs. (33) then imply that the leading corrections to the energies E^ for contact interactions 
are linear in the effective range Ve- 

We have also studied the e-dependence in the three-body sector. In the left panel of 
Fig. 6, we show the corrections from finite e for finite A^. We plot the deviation of results 
for smeared interactions with e 7^ from contact-interactions results with A^ = 70. We pick 
out the unitary limit for contact interactions, which corresponds to E^'^^ = 0.5. Please note, 
that we keep the two-body ground state constant while e is changed. For finite e, the system 
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is not exactly in the unitary limit anymore. For small e there is a linear dependence in the 
double-logarithmic plot. The fit yields a slope round about s = 2 for the state presented, 
which implies a quadratic dependency of the corrections in e for e — )■ in agreement with 
the behavior in the 2-body sector. 




100 



E 



Figure 6. Left panel: Corrections in energy eigenvalues as a function of e with A = 2>, E^'^^ = 0.5, 
E^^^ = 0.5. Right panel: Corrections in (E'^^) — E'^^'^)/ E^'^'^ for the unitary limit with A = 3. 

Finally, there are corrections from a mismatch in renormalization energies E^'^^ and E^^^ . 
In an effective theory, only states with excitation energies small compared to the cutoff 
scale of the effective theory can be described. Clearly, the renormalization energies must be 
chosen in this energy range. For E^"^^ ^ E^^\ we expect errors governed by {E'-^^ - E^^^) 
from the mismatch in the renormalization energies. For contact interactions (e = 0), we have 
investigated the errors in {E^'^'^ — E^^'>) numerically. For this purpose, we chose the unitary 
limit and varied E^^\ The right panel of Fig. 6 shows the deviation of the exact results from 
results in the model space for various {E^'^^ — E^^'^). As expected, the corrections grow with 
increasing \E^'^^ — E^^^\. Thus in practice, it is desirable to choose E^'^'^ ^ E^^^ to minimize 
this type of errors. 



III. RESULTS 

We will now use these insights to study the energy spectra of A-boson systems for A = 
3, 4, 5, 6 in a trap. We will follow two strategies to obtain the spectra from numerical 
calculations: 



1. use contact interactions and extrapolate in 1/ ^/N + {A — 1)3/2 as in our previous 
work [26], and 

2. use smeared contact interactions and extrapolate the converged results in linearly 
in the width parameter e. 

The difference between the two methods will be used to estimate the errors in our calculation. 
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A. Three identical bosons 



At first, we consider the 3-boson sector. In the left panel of Fig. 7, the eigenvalues of 
the first excited L'" = state are shown as a function of the cutoff parameter iV < 70 for 
various smearing parameters e in oscillator lengths /3. The corresponding renormalization 
energies are E^"^^ = 0.5 and E^^'^ = —1. For small model spaces the results are identical 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 2 ^ — — ' • • • • ^ 

^,,^^3,1/2 0.1 0.2 0.3 0.4 0.5 0.6 

MP] 



Figure 7. Left panel: Eigenvalue for the first excited L'^ = 0^ state with A = 3 as a function of N 
with ^(2) = 0.5 and S^^) = _i foj. 

various e. Right panel: The converged eigenvalues as a function 
of e. L'^ = , E^"^^ = 0.5 and E^^'^ = — 1. The points at e = denote the exact eigenvalues from 
[16]. The dashed lines are linear fits to the values for e < 0.3 used to extrapolate to e = 0. 

for all interaction widths. With increasing interaction width e the eigenvalues converge for 
smaller values of N. The converged eigenvalues do not coincide with the results for contact 
interactions. However, as discussed in the previous section, the parameter e can now be used 
as a extrapolation parameter instead of A^. 

The converged eigenvalues are given as a function of e in the right panel of Fig. 7 for the 
5 lowest L'^ = excited states. The solid lines are linear extrapolations of the data points 
for e < 0.3 and the results agree with the exact values known for contact interactions in 
the unitary limit [16] within 3% errors. The error bars of the data points are estimated by 
the difference of the eigenvalues related to the two highest cutoff parameters A^, i.e. with 
these two eigenvalues the exact one is approximated by linear extrapolation and the error 
is estimated as the difference between the state with the highest cutoff parameter and the 
estimated one. In Fig. 7 the two model space sizes correspond to iV = 70 and N = 68. 

In Fig. 8, we show the spectra of Efimov-like states in the unitary limit (E^"^^ = 0.5) 
for E^^^ = 0.5 and E^^^ = —4 in order to illustrate the significance of the mismatch in 
renormalization scales. With a linear extrapolation the exact values for E^^^ = —4 are 
undershot systematically. For larger E^^^ the extrapolation in achieves better results 
than a linear extrapolation in e. Our method clearly has larger errors for E^^'' = —4. 

Finally, we give in Fig. 9 a non-unitary example with E^"^^ = —2. In the left panel the 
eigenvalues for the renormalization energy = -1 are shown as a function of the width 
e and in the right panel for E^^'^ = —5. The results are linearly extrapolated like in the 
unitary limit. At e = the results for contact interactions extrapolated with a quadratic 
polynomial in 1/ y/N + 3 are added. The uncertainties from the extrapolation are estimated 
conservatively as the energy shift from the last calculated eigenvalues to the extrapolated one. 
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e [|3] E m 

Figure 8. Spectrum of Efimov-like L'^ = 0"*" states with ^ = 3 in the unitary limit as a function 
of e. At e = exact eigenenergies are drawn in. Left panel: Renormalization energy ^(3) = 0.5. 
Right panel: Renormalization energy = -4. 

Inside relative errors of 3.5% referred to the ground state with eigenenergy —5 respectively 
— 1 both results coincide. Our method clearly has larger errors for E^^^ = —5. 




E [PI E m 



Figure 9. L'^ = 0^ eigenvalues for j4 = 3 as a function of e with E^"^^ = —2. The points at e = 
denote the extrapolated results for contact interactions. The lines are linear extrapolations to e = 0. 
Left panel: E'-^'> = -1 Right panel: E^^^) = -5 



B. Four identical bosons 

In Ref. [26], we published results for a 4-boson system with contact interactions. Here 
we revisit these calculations using smeared interactions. For instance, calculations for the 
unitary limit E'^'^^ = 0.5 with E^^'^ = —1 with contact interactions e = are compared to 
calculations with e 7^ 0. We can reach cutoff values up to = 26 which is significantly 
smaller than in the three-body sector. In Fig. 10, the eigenenergies for the ground and 
for the first excited states are shown for various model-space sizes. The solid lines are 
extrapolations with quadratic polynomials in 1/ a/A^ + 9/2. As in the three-body sector the 
eigenvalues converge in A^ earlier for larger smearing parameters e. However, for smeared 



15 



-3.5 n . , , . , r- 

\ 

-4.5 

If 

lU 

-5 - 

□ 

□ 

□ 

-5.5 ^ xxxxBxx X X X 

□ * 

-6 L i ' ' — ^ 

0.05 0.1 0.15 0.2 0.25 0.; 

1/(N+4.5)^'^ 

Figure 10. Four-boson L'^ = 0+ states with E^'^'^ = 0.5 and E'^'^'> = — 1 for various smearing 
parameters e as a function of I/^JN + 4.5. Left panel: Ground state. Right panel: First excited 
state. Solid lines are polynomial extrapolation. 

interactions the eigenvalues of the ground state increase at first until they reach a maximum 
and begin to decrease afterwards. It is conceivable that for contact interactions and e = 
0.1 the model spaces are too small in order to see the decrease of the eigenvalues. Thus, 
our extrapolation given in the left panel of Fig. 10 could have large systematic errors. 
Since contact interactions can behave quite differently from smeared interactions, however, 
a definite answer requires calculations at larger cutoffs which are beyond the scope of this 
work. If the renormalization energies E^'^'^ and E'^^'> are chosen equal, the non-monotonous 
behavior in e does not appear. It becomes more severe as the mismatch between E^'^'^ and 
E'^^'> is increased. The excited states do not show the non-monotonous behavior in e at all. 
Their eigenvalues decrease monotonously with e and for large e they converge inside the 
model-space sizes considered here. 

In Fig. 11 (left panel) we show the converged or extrapolated eigenvalues as a function 
of e for the three lowest states with renormalization energies = 0.5 and E^^) = -1. 
In comparison to the three-body sector, we can diagonalize the Hamiltonian only in small 
model spaces (up to = 26). Thus, just until e = 0.35 the eigenvalues have small error bars 
and we extrapolated the states linearly to e = 0. At e = the eigenstates extrapolated via 
A^ with a quadratic polynomial in 1/ \J N + 9/2 and their estimated uncertainties are shown. 
For e < 0.25 the eigenvalues for the ground states cannot be determined with acceptable 
accuracy. The extrapolated results differ significantly and are at odds with each other for 
the ground state. This is due to the non-monotonic behavior of the ground state energy 
with e discussed above. For the excited states this problem is not present and the results 
of the two extrapolations are compatible. In the right panel of Fig. 11, we show the results 
for = 0.5. bmce 

^(2) = ^(3) 

one expects less uncertainties for a fixed e than for 
the case -E*^^-* 7^ E^'^^ (see section II D). In particular, the ground state energies show no 
non-monotonous dependence on e and the extrapolations in e and A^ are consistent. 

C. 5 and 6 identical bosons 

We now apply our combined extrapolation technique to perform exploratory calculations 
of the spectra of 5 and 6 identical bosons. In order to keep the uncertainties as small as 
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e m t [(5] 

Figure 11. Lowest three energy states L"" = 0+ for ^(2) ^ 

as a function of e for A = A. Left 
panel: E^^^ = —1. Right panel: E^^^ = 0.5. The eigenvalues are linearly extrapolated. At e = 
extrapolated eigenvalues for contact interactions drawn in. 



possible, the renormalisations energies are chosen to be identical, i.e. E^'^'^ = E^^'^ = 0.5. In 
particular, this removes the problem of the non-monotonous behavior of the energies as a 
function of e discussed above for the four-body system. The cutoff parameter is = 20 for 
5 bosons and = 16 for 6 bosons. 

Our results for the three lowest energy L'^ = 0^ states with different smearing parameters 
e are shown in Fig. 12. Due to the small model spaces, the uncertainties are significantly 
larger than for three and four bosons. The eigenenergies for contact interactions are extrap- 
olated in with a polynomial of second order in 1/ N + {A — l)/2. They are shown at 
e = with conservative error bands. These approximated eigenvalues are consistent with 
extrapolations in e for smeared contact interactions inside the error bands. Thus we conclude 
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Figure 12. Lowest three energy states L'^ = 0+ for S^^) = q.S, E^^^ = 0.5 as a function of e. Left 
panel: A = 5 {N = 20). Right panel: ^ = 6 (A^ = 16). The solid lines are linear extrapolations in 
e. At e = the eigenvalues extrapolated with a polynomial of second order in l/y^A^ + {A — l)/2 
for contact interactions drawn in. 



that the combined extrapolation in A^ and e makes calculations in 5- and 6-boson systems 
with moderate computational resources possible. 

In Table I, we have collected the extracted energies of the lowest three Efimov-like states 
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for the unitary limit {E^"^^ = 0.5) and a three-body ground state energy fixed to E^^^ = 0.5 
(of. Figs. 8, 11, and 12). The numbers in the first two columns are the exact values for A = 2 



A = 2 


A = 3 


A = 4 


A = 5 


A = 6 


0.5 


0.5 


-0.1(2) 


-0.9(2) 


-2.3(3) 


2.5 


2.9 


2.7(3) 


2.2(3) 


1.1(2) 


4.5 


5.1 


4.6(5) 


4.3(5) 


3.7(3) 



Table I. Energies of the three lowest Efimov-Iike states in systems with A = 3,4,5,6 for the 
rcnormahzation energies E^"^^ = 0.5 corresponding to the unitary limit and E^^^ = 0.5. The column 
labeled A = 2 contains the three lowest two-body states. All energies are in units of Huj. 

and ^4 = 3 rounded to two digits of accuracy. The values for the columns labeled A — A, 
A = 5, and A = 6 are extracted from our calculations. Here the number in parentheses 
gives the difi'erence in the last digit between the extrapolation in N and e which can serve 
as an estimate of the numerical error. The spectra are much more compressed than the 
corresponding free space results [4-9]. However, a similar pattern can be observed: starting 
from ^4 = 3, every A-hody state is accomponied by a corresponding A -\- 1-body state. A 
second A + 1-body state attached to the A-body state is not present in our calculation. It 
would be interesting to study the spectra also away from the unitary limit for fixed E^"^^ to 
investigate their systematics. However, due to the additional errors from the mismatch in 
E^"^^ and E^^^ this is numerically more challenging. 

IV. CONCLUSION AND OUTLOOK 

In this paper, we have extended our previous work on universal few-body physics in 
harmonic confinement [26]. We investigated A-body systems with up to six particles using 
genuine and smeared contact interactions at leading order in the large scattering length. 

The effective theory is defined within a finite model space with a cutoff on the basis 
functions. Wc perform a detailed study of the running of the effective two- and three-body 
coupling constants g^"^^ and g'^^'> with for contact interactions. In the two-body sector, the 
results of Stetcu et al. [28] for the running of gf*^^-* are confirmed. For g'^^\ the limit cycle 
behavior expected from free space is observed and one complete cycle is traced numerically. 

Moreover, we have performed a detailed analysis of the leading effective theory errors 
using Lepage plots. Our cutoff A^ on the size of the model space implies both an infrared 
cutoff of order N + {A - 1)3/2 and an uhraviolet cutoff of order l/^jN + {A - 1)3/2 in 
the effective theory. The errors in the eigenenergies from these cutoffs for genuine contact 
interactions scale as 1/ -^/N -\- {A — 1)3/2 in agreement with the expectations from the con- 
tinuum case [3]. For smeared contact interactions, we find an exponential dependence on 
\fN (cf. Eq. (28)) in agreement with Ref. [33]. 

An additional error is introduced if there is a mismatch in renormalization energies 
and E^^^ . We investigate these errors for contact interactions numerically and find that they 
are governed by the difference [E^'^'> — E^^^). To improve convergence E'^^^ and E'^^^ should 
be chosen as close as possible. 

Smearing the contact interactions over distances of order e improves the convergence 
with A^ considerably. Obviously, this smearing changes the interaction. In order to extract 
results for genuine contact interactions from calculations with finite e, extrapolations in e 
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are required. These extrapolations show a quaUtatively different behavior for finite N and 
infinite N since the limits N ^ oo and e — > do not commute. 

Based on our findings, the optimal strategy to extract the universal energies for contact 
interactions from numerical calculations in a finite model space when convergence can not 
be reached is to (i) use genuine contact interactions and extrapolate to = oo and (ii) use 
converged energies for smeared contact interactions and extrapolate linearly to e = 0. The 
difference between the two methods can be used to estimate the errors in the calculation. 

We use this strategy to calculate the spectra for A-boson systems with up to six bodies 
and study the case E^'^^ — E'^^^ — 0.5 in detail. We find that starting from A — 3, every A- 
body state is accomponied by a corresponding A + 1-body state. A second A + 1-body state 
is not observed. The combined extrapolation in and e makes the calculations possible 
with moderate computational resources and provides a first step towards calculating trapped 
many-body systems with resonant interactions. On the experimental side, the deterministic 
preparation of tunable few-fermion systems of up to ten particles in a microtrap was recently 
achieved by Serwane et al. [40]. Similar experimental studies of few-boson systems appear 
also promising. 

If there is a mismatch in E^'^^ and E^^\ the ground state energies for four and more 
particles show a non-monotonous dependence on N such that reliable results can only be 
obtained in large model spaces. For the excited states this problem is not present. In the 
future, it would be interesting to get a better analytical understanding of these corrections 
and possibly develop an renormalization group improved matching procedure for this case. 

Finally, it would also be interesting to remove the trapping potential from the Hamiltonian 
(2) and apply our method to universal bound states in free space. For nuclear systems 
in the framework of the no-core shell model such a strategy was applied by Stetcu and 
collaborators [23]. 
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Appendix A: Construction of the Hamiltonian Matrix 

In this appendix, we discuss the construction of the Hamiltonian matrix for spinless 
bosonic wave functions, i.e. symmetric functions under transposition of the bodies. A sym- 
metric basis for the model space of A bodies is constructed iteratively starting from the 
symmetric basis for two bodies. 

a. Symmetric basis 

To construct unnormalised symmetric states we define P*^^-* = (1 + T12) with the transpo- 
sition Tjfe between body i and k. With the oscillator states and the multi- index u we 
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thus have for 2 particles 

[2(pni,n{si) for < even, 
The orthonormal basis for 2 particles with multi-index a thus is 

Wa{si)^^nlm{si) for Z odd. (A2) 

From these states coupled mixed-symmetric 3-body states are constructed as [0''^^Q,®0nz]^- 

Since [H, L] = 0, the matrix elements are independent of M and this quantum number is 
suppressed. 

Accordingly, on the symmetric, unnormalised A-particle states we have 

A- 2 

piA) ^ (1 + j.^^ + . . . + T^a-i)a)P^^-'^ = (1 + E T^A-i)iT^A-i)ATi^A-i) + T^a-i)a)P^^-'^ . 

i=l 

Furthermore, we introduce the transformation matrix 

Ci^rn'^^^W, n',, l'„ L'] = { [W^^, ® J I [^,] > , (A3) 

with multi-indices a' and z from the mixed-symmetric 74-body states <8) ^^n'^v^ to 

the symmetric A-body states ^^^^^j, given by 

^sym "^''"^"^[^j ,n2,l2, L'] 

= {[W^a' ® Mlif\P^^^\[¥^a ® 0n./.]'> , (A4) 
= {[W^a' ® Ml'fli^ + - mA-l)A)\[W^a ® <Pn,l,]') , (A5) 

where z = {a, 77-2, h, L}. 

Our goal is to construct an orthonormal basis of symmetric 74-body states (p^^^^^- By 
diagonalization of the symmetric norm matrix 

J\f[z';z]^{W^,\W^,,), (A6) 

the diagonal matrix D of non-negative eigenvalues is found: 

D = 0-X-0^. (A7) 

Vanishing eigenvalues do not correspond to symmetric states and respective eigenstates 
are eliminated. Thus, we find the non-quadratic transformation matrix i^gym from mixed- 
symmetric A-body states to orthonormal symmetric states: 

B^^-'^-'^iu^; a', n'„ l',L'] = -j^=0\j; z]Aif-'^^^[z; n'„ l',L'] , (A8) 

where a sum over j and z is implied. Accordingly, the orthonormal basis of symmetric states 
is given by 

^iym '^"^[^; ^2, k, L] I [W^^ ® . (A9) 
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b. Explicit calculation of ciym ^^^^ 

The explicit calculation of Ci^^''^"^ is based on the Talmi-Moshinsky-transformation 
and on a change in the coupling scheme. With the operator Eq. (A5) yields 

(AlO) 

In order to calculate matrix elements of the operator T(a-i)a, a transformation to coordinates 
A and jl is performed since T(^a-i)a is diagonal in A and jl (see Fig. 13): 




A-2 
2{A-l) 



2(A-l) 



2(A-1 ) 

A-2 
2(A-1) 



S(A-2) 
S(A-1) 



(All) 



For this purpose the coupling scheme is changed by means of Wigner's 6j symbols [41] and 
in the new coupling scheme the Talmi-Moshinsky-transformation is exploited. 

With the Brody-Moshinsky-brackets {nxlx,nJn; Li2\nili,n2l2) ^_ we find for the wave 

functions in terms of the coordinates A and /2 



(A-2) 



E 



^ni,/i(S(A-2))]^" ® 0n2«2(^(A-l) 

J2 (-1)^-+'^+'^+^ V(2Li2 + 1) (2L, + 1) I 



I2 L L12 



nxlx,n^l^; Li2\nili,n2l2) ^ 4>'^^-'^\'' (S) [(pn^^i^W ® 4>n^i^{f^)\ 



(A12) 



c. Contributions of the 2- and 3- body contact interactions 

In addition to the matrix Cgym, also the 2-body interactions can be calculated in the 
coordinates A and fl introduced in the last section (see Fig. 13). Note that the contributions 
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of all 2-body interactions are equal since the bodies are identical. The 2-body interaction 
between body {A — 1) and A depends only on the distance between these, which is propor- 
tional to A and therefore the interaction is proportional to 5*^^^ (A). With combinatorics the 
number of interactions are determined: There are (^) = ^^^^"""^ two-body interactions. 

Using the recoupling and the Talmi-Moshinsky-transformation one more time the contri- 
butions of all (^) = three-body interaction are determined in a similar fashion: 

For states expressed in terms of the coordinates z7, k and A the interactions can be readily 
calculated (see Fig. 14). The corresponding transformation is 




S(A-3) 



(A13) 




(A13) 



A N-1 



N-3 N-2 




N A N-1 



Figure 14. Talmi transformation from coordinates S(A-3) ^-^d fl to coordinates z7 and k 

Appendix B: Smeared contact interaction 

In this section we give the derivation for the matrix elements 



1 , /I - eS"! , /I - e^\< 



{nih\VG\n[l[) = ^^^^^<PnAO)[^) 0n;o(O)(^) ^5.,o5r,o • (Bl) 



For the interaction 



we have 



\Vg\s) 



1 _ _ s'^ 

e 2? e 



(27re2) 



(B2) 



1 /' I' s'^ 



(B3) 



Moreover, 



_2 /"OO 2 1"^ /Is 

dh 0„z(s)e-^ = 74^ / dr r^e-^R„o{r) = N^o / dr r^e~^^^'^^'>''' rn^\r^) , (B4) 

Jo Jo 
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with the associated Laguerre polynomials and the normalization A^„o = '\J'^^^^^u,- 

By the substitution of = x and with — |(1 + we get with the expansion of the 
Laguerre polynomials in monomials 

47rAr,o (-1)' r(n+|) 1 



2 i\ r(n-i + l)p3 



4vriV„or(n+i)l^,.> ^^^^ 



i=0 ^ 



2"+2^!p2n+3 

and replacing p and A^„0) we find 



4.iV.o(2n + l)!!^^,_^^„^^ (B8) 



/ 1 /(2ni + l)!! /(2n; + l)!! /l-e2\"Vl-e'\"^ ^ 



= (rTi^'^-o(0)(^)"Vnlo(0)(^)"'^5.,o5^^^ . (BIO) 

Appendix C: Effective range expansion for smeared contact interactions 

For separable potentials 

(sills') = ^(2)^(s>(s-'), (CI) 

we get for the T-matrix 

i,^T{z)\^) = - B{z))^ \{p)v{0) with B{z) := J (C2) 

from the Lippmann-Schwinger equation 

T{z)^V + V{z-Ho)-'T{z) . (C3) 
Here, v{p) denotes the Fourier transform of a;(s) 

v{p)^ J dh e'^uj{s) . (C4) 

For separable central potentials there is only s-wave scattering and we find for the scat- 
tering amplitude fo{p) with the effective range expansion 

fo{p) = 1.1 2 ! /o^ 4^ ~ 

= -i- hm (p-lT(pV2 + zA) \^) (C6) 
= [-2n\v(p)\-'{l/g^'^ - B(py2 + iA))]-' . (C7) 
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For the smeared contact interaction we have 



uj{s) 



(27re2) 

-j^ poo 



3 ^ ) 
2 



v{p) = e 2 , 



dk 



Vke-'^'' 



and the imaginary part of the denominator in /o(p) can be calculated with 



J a 



dx 



X 



It follows for the imaginary part 



fo 27r^ 



and for the real part 



5^1 
fo 



27r 



1 

TT 



^ — c 



dq 



g,2 _ p2 



A(p2) 



Moreover, it can be shown, that 

where the Dawson function is defined by 



2epF(ep) 



F{x) 



f dye^'^x- '^x^ + 0{x^) . 
Jo 3 



(C8) 



(C9) 



(CIO) 



(Cll) 



(C12) 



(C13) 



The function A{p'^) as well e^^^^ is expanded in p"^. We will be interested in the scattering 
length a and the effective range r, thus the expansion is terminated at 0{p'^): 



- + \rp^ + 



2n 



1 vTT / / 2 o o 

--^ l-2ep[ep ep' 

TT e V ^ 3 



1 + eV 



0{p') , 



_ (2ti 1 . 1 / 2e 



(CM) 
(C15) 



Note, that the momentum p is the Jacobi momentum. The relation between the oscillator 
length for canonical relative coordinates and Jacobi coordinates is just a factor of -\/2 in the 
two-body sector and the oscillator lengths differ with the same factor. Therefore all lengths 
must be rescaled with -\/2 in order to get the canonical effective range expansion. 
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